arXiv:1501.04065v2 [hep-th] 2 Mar 2015 


Numerical determination of OPE coefficients in the 3D Ising model from off-critical correlators. 

M. Caselle\ G. Costagliola \ N. Magnoli^ 

^ Dipartimento di Fisica, Universitd di Torino and INFN, Via P. Giuria I, 10125, Torino, Italy. 

^ Dipartimento di Fisica, Universitd di Genova and INFN, Via Dodecaneso 33, 16146, Genova, Italy. 

(Dated: March 3, 2015) 

We propose a general method for the numerical evaluation of OPE coefficients in three dimensional 
Conformal Field Theories based on the study of the conformal perturbation of two point functions in the 
vicinity of the critical point. We test our proposal in the three dimensional Ising Model, looking at the 
magnetic perturbation of the < CT(r)(j(0) >, < (T(r)e(O) > and < e(r)e(0) > correlators from which we 
extract the values of CJe = 1.07(3) and CIj = 1.45(30). Our estimate for C^e agrees with those recently 
obtained using conformal bootstrap methods, while Cl^, as far as we know, is new and could be used to 
further constrain conformal bootstrap analyses of the 3d Ising universality class. 


INTRODUCTION 


In the past few years a renewed interest in the Confor¬ 
mal Bootstrap approach IT]] to three dimensional Confor¬ 
mal Field Theories (CFT) led to a set of remarkable re¬ 
sults iH-iIii on the universality classes of several 3d sta¬ 
tistical models. Among these models a particular attention 
has been devoted to the 3d Ising case, both for its phys¬ 
ical relevance and for the fact that in this case very pre¬ 
cise numerical estimates exist for the scaling dimensions 
which allowed highly non trivial tests of the confor¬ 
mal bootstrap results. It would be interesting to perform a 
similar comparison also for the Operator Product Expan¬ 
sion (OPE) coefficients 112 ih . but for these constants Monte 
Carlo estimates with an accuracy comparable with that of 
the scaling dimensions are still lacking. As an attempt to 
fill this gap we propose in this paper a general strategy for 
the numerical estimate of the structure constants and, as 
a proof of concept, we evaluate and in the Ising 
case. Our estimates should be considered mainly as an ex¬ 
ploratory study on the effectiveness of the procedure with 
a limited use of computer power. We plan in a forthcoming 
paper to obtain more precise estimates for these and other 
OPE coefficients. 


It is well known that estimating structure constants is 
much more difficult than estimating scaling dimensions. 
The most naive approach would be to extract them directly 
from three point functions (i.e. the connected correlators 
of three operators) at criticality. However it is easy to see 
that the subtraction in the connected correlator of the two 
point components leads to a mix of contributions with sim¬ 
ilar exponents which makes it very difficult to extract the 
sought for OPE coefficients. In this letter we propose an 
alternative method based on conformal perturbation which 
allows to overcome this problem. The main idea is to study 
the perturbed two point functions using the OPE to rewrite 
them as sums of perturbed one point functions and then use 
the known scaling behaviour of one point functions and 
structure constants to write the perturbed two point func¬ 
tion as an expansion in powers of the scaling variable. Eol- 
lowing 01 it is always possible to write the short distance 


behaviour of two-point correlators in the vicinity of a criti¬ 
cal point as a short distance expansion: 


< O,(r)Oj(0) >h= < 0^(0) >,, (1) 

k 


where < .. denotes the expectation values with re¬ 
spect to the perturbed action, the Oi represent a complete 
set of operators of the conformal theory and the (r, h) 
are the Wilson coefficients. While the one point functions 
< Ofc(O) cannot be determined from the knowledge of 
the critical correlators and take care of the non-perturbative 
physics of the model, the Wilson coefficients can be ob¬ 
tained from integrals of critical correlators. Moreover it 
was shown in ll22ll that the order derivatives of the Wil¬ 
son coefficients with respect to the perturbing parameter of 
the theory are always infrared finite. The structure con¬ 
stants can then be easily extracted from the coefficients of 
the short distance expansion (see 12214241 for further de¬ 
tails). In practice one is limited for numerical reasons to the 
first few coefficients (typically the first two or three) which 
however (combining all the possible two points functions 
and relevant perturbations) are enough to obtain all the 
structure constants among relevant operators. 

As it is easy to see the method is based on very gen¬ 
eral features of CETs and OPEs. It can be applied to 
any d-dimensional CET and in fact it was applied a few 
years ago to the 2d Ising model perturbed by a magnetic 
field 11^ l24 ], to various perturbations of the Tricritical 
Ising model and then used to identify the signatures of 
subleading irrelevant operators in the perturbed two point 
functions J^. In all these cases the structure constants 
were already known from the exact solution of Belavin 
Polyakov and Zamolodchikov 0 and the main emphasis 
was on comparing the perturbative results with the numer¬ 
ical simulations. 


The aim of this paper is to apply this strategy to the three 
dimensional case. In particular we shall address, as an ex¬ 
ample, the magnetic perturbation of the 3d Ising model. 
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ISING MODEL PERTURBED BY A MAGNETIC EIELD h 

The action of the continuum Ising model in 3D per¬ 
turbed hy a magnetic field is: 


where < i,j > indicates nearest neighbors sites in the 
lattice which we assume to be a three-dimensional cu¬ 
bic lattice of size L. Fixing /3 at its critical value = 
0.2216544 |28| and defining hi = f3cH, we have: 


S = Scft + h j a{r)dr. (2) 

where S^ft is the action of the conformal field theory which 
describes the model at the critical point. 

At the conformal point there are two quasi-primary rele¬ 
vant fields cj and e, whose dimensions are given by A^- = 
0.51814(5), A, = 1.41275(25) Q. 

The VEVs acquire a dependence on h, that can be fixed 
using renormalization group arguments: 

< o > = A„h'^ (3) 

< e > = (4) 

In order to simplify notations let us introduce the scaling 
variable t = and an adimensional version of the 

structure constants and their derivatives: 

C’^. = lim Q) 

r->oo ^ 

due'll = lim 

With these definitions the short distance expansion of a and 
e correlators becomes: 

< a(r)a(0) > = Cl^ + (5) 

+ A^d,C:j^+^ + Oie) 

^ + dnClt (6) 

+ + 0(f^+") 


Z = '^ (9) 

cri=±l 

It is natural to define the lattice discretization of a as 
(T* = -p ^ ■ cJi, so that its mean value coincides with the 
magnetization. In a similar way one can define the energy 
operator e* as e* = Aj J2<ij> ~ where ecr is 
the energy at the critical point, so that the mean value of 
coincides with the singular part of the internal energy. 

It is important to notice that measuring the mean values 
(I3I4I) and the two point functions (I5I6I7II on the lattice we 
find the lattice versions of the amplitudes A and structure 
constants C (which we shall denote in the following with 
the index 1: A’’ and C*). To relate these constants with the 
continuum ones we must first fix the relative normalization 
of (j*, and hi with respect to a, e and h. The simplest 
way to do this is by comparing the correlators at the critical 
point. In fact from 

< CJiCTj >= I 

we get (assuming the standard normalization of the two 
point function in the continuum) cj* = By consid¬ 

ering the e* correlator we get similarly e* = R^e. Finally, 
hi is fixed by matching the lattice and continuum perturba¬ 
tion terms: hi = R~^h. 

Using R^ and R^ we can easily obtain continuum ampli¬ 
tudes and OPE constants from the lattice ones. Eor example 

from the equation dU) we get = R^'^ ReA^ and from 
this and the equations (l5]l and (fTOl l we have the rescaling 
for the structure constant: CijR. = Ct„. 

(T(T t (7(7 


< e(r)e(0) > = Cl + A,Clt^ (7) 

+ A,di,Clt^^^ + o{e) 

where A/i = 3 —Ag. = 2.48186(5) is the dimension of the 
magnetic field and with the standard GET normalizations 
we may set Cl = = 1. 

Conversion from lattice to continuum normalizations. 
The partition function of the 3d Ising model is given by: 

cri=±l 


MONTE CARLO SIMULATIONS 
Estimating and Re 

As a first step we extracted the values of R^ and Re from 
a finite size scaling analysis of the two point correlators 
at the critical point. The table below summarizes our re¬ 
sults. We are not aware of independent estimates of these 
quantities except for ecr which was estimated in lE^ as 
(cr = 0.3302022(5) in good agreement with our result. 
The quoted uncertainties combine both statistical and sys¬ 
tematic errors. It should be noticed, in view of possible 
future improvements of our analysis, that the main source 
of error was due to the infinite volume extrapolation of our 
results (we used lattices in the range 120 < L < 200). 
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Table I. Rescaling and other useful constants 


Rcr 

0.550 (4) 

Re 

0.237 (3) 

Ai 

1.0125 (5) 

Ai 

0.608 (1) 

^cr 

0.330213 (12) 


Perturbed correlators 


Table II. Results for the structure constant C%a obtained from 
the spin-spin correlator. We report systematic errors in square 
brackets. 


hi 

^min 

max 


2 • 10“'" 

6 

20 

1.06 (3)[7] 

10“® 

5 

35 

1.04 (2) [7] 

0.85 • 10"'^ 

5 

35 

1.06 (3)[8] 

0.75 • 10"® 

5 

40 

1.07 (3)[8] 


To extract the structure constants we then estimated the 
same correlators in presence of a small magnetic field hi. 
The choice of hi is constrained hy two main requirements: 
hi should he small enough to keep the correlation length 
^ as large as possible, at the same time it should not he 
too small to avoid finite size effects. The optimal choice is 
thus fixed hy the maximal lattice size which we could sim¬ 
ulate with the computer resources at our disposal, which 
was fixed to he L = 200. Given this constraint the optimal 
range turned out to he 2 • 10“^ < hi < 0.75 • 10“^. In 
this range the correlation length spans from about 20 to 40 
lattice spacings, allowing to sample a sufficient number of 
different distances of the correlators. 

For all values of h we performed our simulations on lat¬ 
tices of size L = 200 using a state of art Monte Carlo 
algorithm with about 10^ configurations for each simula¬ 
tion. We fitted the correlators with eq.s ([5] |6] I?]), setting 
= Cl^ = 1, using A^., Ag as fixed inputs and 

keeping as only free parameters the structure constants. It 
turned out that in all correlators the last terms quoted in 
eq.s (|5l0 and the last two terms in eq. |6]were negligible 
within the errors and in all three cases we ended up with a 
linear fit with only one free parameter. A major source of 
uncertainty in our estimates is the systematic error due to 
the uncertainty in the estimates of and We quote 
separately in the following these errors (which we report in 
square brackets) from the statistical ones. Due to the dif¬ 
ferent scaling behaviour, the relative size of statistical and 
systematic errors is different for the three correlators. For 
estimated from < (T(r)cj(O) > the two errors are of 
similar magnitude. Cl^ extracted from < e(r)e(0) > is 
dominated by the systematic error while for extracted 
from < (T(r)e(O) > the systematic error is negligible (due 
to the fact that the critical correlator vanishes in this case). 

The tables below report our results for the structure con¬ 
stants. Tmin and rmax dcnotc the range of distances in¬ 
cluded in the fit. In all cases we chose = 5. We 
verified that in all three cases this was enough to elimi¬ 
nate lattice artifacts and to neglect short distance sublead¬ 
ing contributions in the correlators. 

Looking at tab.s (IIIIIVI) we see that, within our statistical 
errors Combining the results listed in tab.s 

(IIIIIIIIIVI) we quote as our final estimate for the two struc¬ 
ture constants = 1.07(3) and = 1.45(30). The 
first value is in good agreement with a recent conformal 


Table III. Results for the structure constant (71^ obtained from the 
energy-energy correlator. We report systematic errors in square 
brackets. 


hi 

min 

^ max 

Cle 

2■10"® 

5 

16 

1.47 (10) [30] 

10"® 

5 

18 

1.42 (10) [50] 

0.85 - 10"® 

5 

18 

1.43 (12) [50] 


bootstrap calculation ifT^ : (C^J^ = 1.10636(9) H^l- We 
are not aware of any estimate for the second constant 
The very fact that it is different from zero is rather non 
trivial. In fact in the two dimensional model = 0 as a 
consequence of dual symmetry. Our calculation shows that 
indeed this is not a generic property of the Ising model but 
is a specific feature of the two dimensional self-dual case. 


CONCLUDING REMARKS 


In this paper we proposed a general strategy for the 
Monte Carlo estimate of OPE coefficients in d-dimensional 
spin models. As a proof of concept of our method we per¬ 
formed an exploratory study in the 3d Ising case and found 
a value of in agreement with the known conformal 
bootstrap results and a preliminary estimate for which 
could be used as input for further constrain existing con¬ 
formal bootstrap calculations. Similarly, the knowledge of 
these structure constants could be used to constrain the lin¬ 
ear response dynamics at finite temperature of Conformal 
Quantum Critical systems JUj]. 

We see a few possible improvements of our approach: 

• The systematic uncertainty which is the major source 
of error in Cl^ could be improved, following the ap¬ 
proach proposed in lE5l . looking at different realiza¬ 
tions of the 3d universality class with improved scal¬ 
ing behaviour. 

• It would be useful to combine the magnetic perturba¬ 
tion that we studied in this paper with other types of 
perturbation. 


It would be interesting to extend the method also to 
boundary CFTs ® and to CFTs containing confor¬ 
mal defects 133. 341 


We plan to address these issues in a forthcoming publica¬ 
tion. 
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Table IV. Results for the structure constant Cje obtained from the 
spin-energy correlator. 


hi 

'^min 

'^max 

^cre 

2•10“^ 

5 

16 

1.10(3) 

10“® 

5 

16 

1.07 (4) 

0.85 • 10“*^ 

5 

14 

1.07 (5) 
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